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Abstract The variational local moment approach (V-LMA), being a modification of 
the method due to Logan et ai, is presented here. The existence of local moments is 
taken from the outset and their values are determined through variational principle 
by minimizing the corresponding ground state energy. Our variational procedure 
allows us to treat both fermi- and non-fermi liquid systems as well as insulators 
without any additional assumptions. It is proved by an explicit construction of the 
corresponding Ward functional that the V-LMA belongs to the class of conserving 
approximations. As an illustration, the V-LMA is used to solve the multi-orbital 
single impurity Anderson model. The method is also applied to solve the dynamical 
mean-field equations for the multi-orbital Hubbard model. In particular, the Mott- 
Hubbard metal-insulator transition is addressed within this approach. 



1 Introduction 

The single impurity Anderson model (SIAM) is one of the most investigated mod- 
els in condensed matter physics HI. This model is regarded as a prototype to un- 
derstand and describe: i) properties of metals with magnetic atoms |T6l, ii) charge 
transport through quantum dots 1241 . iii) Mott-Hubbard metal-insulator transitions 
(MIT) within the dynamical mean-field theory (DMFT) l23ll6l[T4ll28l[30ll27l. and 
iv) a crossover between weak and strong coupling limits and confinement phenom- 
ena. The SIAM consists of a term describing band electrons coupled by hybridiza- 
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tion to a term corresponding to a single impurity where the local Coulomb inter- 
action is taken into account HI. In the featureless hybridization limit the SIAM 
is solved exactly within the Bethe ansatz or conformal field theory techniques so 
the ground state and the whole excitation spectrum as well as thermodynamics are 
exactly known lfT6l . Unfortunately, these methods cannot in practice provide dy- 
namical quantities, for example one-particle spectral functions or dynamical sus- 
ceptibilities, for all interesting energies. Also the (asymptotic) exact solvability is 
not possible for a general hybridization term. 

For practical applications of the SIAM one has to rely on either a numerically 
exact or an analytical but approximate solution. Numerically exact methods, like 
the numerical renormalization group (NRG) IS) or the determinant quantum Monte 
Carlo (QMC) 1 14| are very time (CPU) consuming. In particular, the CPU is very 
long when the number of orbitals is large in the NRG case and when the temperature 
is low in the QMC case. Also to extract dynamical quantities is a rather tricky task 
ifm . Reliable analytical methods are therefore needed. One of such methods, which 
recovers properly both weak and strong coupling limits, is a local moment approach 
(LMA) invented recently by Logan et al. 1261 . 

The LMA is a perturbative method around an unrestricted Hartree-Fock solution 
with broken symmetry, i.e. with a non-zero local magnetic moment. The broken 
symmetry is restored at the end by taking the average of the solutions corresponding 
to different directions of the local magnetic moment 1261 . 

In the present contribution we describe the LMA method and our implementation 
of it, which is different from the original one ll26l by the way of how the value of 
the local moment is determined. Namely, we use the variational principle demanding 
that the ground state energy is minimized by the physical value of the local moment. 
Therefore we use the name variational local moment approach (V-LMA) for this 
method. Such a procedure allows us to easily generalize the V-LMA for multi-orbital 
models as well as for finite temperatures and systems with disorder IS] |7] |9] (TO) . We 
also discuss the Luttinger-Ward generating functional for the V-LMA and claim that 
this method belongs to the class of conserving approximations. The application of 
LMA for studying the electron flow through quantum dots and the Mott-Hubbard 
MIT is addressed at the end of the contribution. 



2 Local moment method in one orbital SIAM 

The single impurity Anderson model is given by the Hamiltonian 



^^SlAM — He + Hi,„p + ///lyfo, 



(1) 



where the conduction electrons are described by 



(2) 
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where is an energy (a dispersion relation) for an electron in a state k and spin a = 
±1/2, the impurity electrons with the local Coulomb interaction U are represented 
by 

Hi,„p = Y,i^d + Und-a)"-da, (3) 
a 

with n^ia = dada, and the hybridization between conduction and impurity electrons 
is 

///n* = LK4cka+/^-C-)- (4) 
k,<T 

All local (on impurity site) properties are expressed by the hybridization function 

and not by et and Vk separately. This can be proved by tracing out the non- 
interacting conducting electrons. 



2.1 Mean field solution of the single impurity Anderson model 

The Hartree-Fock mean-field solution of the SIAM is obtained by factorizing the 
interacting term n^/i/vj. ~ {nd\)nd[ +«f/t(«rfi) — {nd\){nd{) CI- For the interaction 
U above Uc and corresponding impurity electron densities «j the mean-field solution 
is unstable toward the local moment formation with non-zero moment /i = («^/|) — 
(njj). The solution is doubly degenerate because of two equivalent directions of 
the local moment /i = ±|Ai|, which give the same energy of the system. The local 
(impurity) Green function within the Hartree-Fock solution is 

where the static Hartree-Fock self-energy — U {iia) and 5 — > 0^. Since there are 
in principle two possible signs of the local moment, there are two different possible 
Hartree - Fock Green functions denoted by G^{co)^^ and G^{co)^^ that differ only 
by the sign of the local moment and depend parametrically on its value |/i |. 

The fundamental deficiency of the Hartree-Fock approximation is that it leads to 
a broken symmetry solution which cannot persist in the thermodynamic limit, i.e. 
a single impurity cannot lead to the magnetic solution in the infinite system. Also 
this solution does not recover the singlet ground state known from the exact Bethe 
ansatz solution. Nevertheless it turns out to be useful as a starting point in the further 
perturbative calculation combined with the symmetry restoration. 
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2.2 Two self-energy description 



The two Hartree-Fock Green functions Ga'^{(o)^'^ are used in the time-dependent 
many-body perturbation expansion. Within the random phase approximation (RPA) 
the polarization diagrams are 

and correspond to spin flip processes as represented by the Feynmann diagrams in 
Fig. [I] 




Fig. 1 The frequency dependent part of self-energy expressed as the RPA series around the broken 
symmetry Hartree-Fock solution. The transverse spin polarization bubbles constitute a geometric 
series which can be summed up to infinity. 



For each type of the mean-field solution G^'^{co)^^ we have the corresponding 
self-energy 

eUco) = eS' + u'I ^^M§t{<o')GU(o~co'r (8) 

depending on frequency and parametrically on as well. The full RPA-Green 
functions Ga'^{(o) are constructed by using the Dyson equation separately for A 
and B solutions. Note that 0^ (ft)) depends parametrically on still unknown |/z |. 



2.3 Symmetry restoration ansatz 

To restore the spin-rotational symmetry Logan et al. Il26l proposed the following 
ansatz for the full symmetrized Green function 

Ga{(0) = \{Gi{(0) + Gl{(0)). (9) 

Within the LMA the physical Green function is an average of the two solutions with 
equal probabilities. Although each G^'^(c(j) is determined within the renormalized 
perturbation scheme the final Green function turns out to capture nontrivial non- 
perturbative physics as was shown by Logan et al. |j26i and is also reproduced below. 
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In particular, the LMA is able to recover the Kondo peak in the spectral function 
correctly with the exponential width. 



2.4 Determining the value of local moment 

The value of the local moment is a free parameter and must still be determined. 
In the original approach, Logan et al. ||26]| imposed the Fermi liquid condition to 
determine at zero temperature. This condition might be too restrictive at finite 
temperatures or in the multi-orbital cases. Therefore we decided to find the physical 
solution to the problem by minimizing the relevant thermodynamic al potential with 
respect to |/J | ||20| . At zero temperature the relevant potential is just the ground state 
energy of the system, i.e. 



where in the case away of half-filling the particle density n must also be determined. 
The variational method reproduces the Fermi liquid properties where they are ex- 
pected. 



2.5 Ground state energy in the Anderson impurity model 

The ground state energy of the SIAM is given by Eq ~ (0|//|0). This quantum- 
mechanical average consists of two parts: the bulk, which is proportional to the sys- 
tem volume and is independent of the local moments, and the impurity part, which 
depends explicitly on \fx\. The impurity part of the ground state energy, expressed 
by the local Green function Ga{(0) and the hybridization function A {(£>), is equal to 



where the contour integral is over the half circle in the upper complex plane. 



2.6 LMA as a conserving approximation 

According to Kadanoff and Baym [S) any approximate theory is conserving if there 
exists a Luttinger-Ward functional 4>[G] for this theory. It is necessary that this func- 
tional: i) is universal, i.e. it dependents only on the full propagator Ga{(0) and not 
on the atomic properties of the system and ii) has a functional derivative with re- 



Sphysicai = minEGiil.n), 



(10) 



EH 




(11) 
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spect to Ga{co) which is by definition equal to the self-energy of the system. It can 
be shown ||20| that the LMA is a conserving approximation and we can construct 
explicitly the Luttinger-Ward functional 

4>[G] = <I>[G^G^] = i {<fipA+<PA) + \Tr\ogGiG', + 

-^rrlogQ(G^, + G^)), (12) 

where Tr = TY,a L/ra,, ^'^"^ the functionals are represented diagrammati- 
cally by the RPA diagrams with G'^^{(o) respectively. The constraint that G = 
5 {Ga + Gb) must be satisfied. Finally, the free energy functional is given by 

Q.[G] = ^[G] + Tr\ogG-TrEG (13) 

and the stationarity condition 5Cl [G]/5G gives the Dyson equation and the physical 
solution for G. 

The fact that the LMA is a conserving approximation, as we proved above, makes 
this theory reliable in describing correlated electron systems, in particular in the 
intermediate regimes of parameters. 



3 Local moment approach for the multi-orbital SIAM 

In reality the magnetic impurities in metals are atoms with partially field d- or f- 
orbitals. Such orbitals have degenerate levels. Even when a particular environment 
which decrease the symmetry and leads to and t2g split levels, partial degeneracy 
between orbitals remains. The appropriate model to describe such situations is the 
multi-orbital single impurity Anderson model. It describes a single impurity with 
many orbital levels a, which can be degenerate or split depending on the single- 
body matrix element £«. In this case the electrons can interact via direct (density- 
density) type of the interaction and via the exchange (Hund) interaction. Microscop- 
ically, the single impurity Anderson model with many orbital levels is given by the 
Hamiltonian: 

^SIAM = 51 (^a + Uana.a)na.a + X! X! ^-^^aa') naan^a' + 

a, a a.a'a^p 

+ L Vka (di^c^^^+clJ„^ + ^ Ek^a^ka' (14) 
k,a,a k.a 

where the direct U and U' as well as exchange J interactions between the electrons 
of spin (7 and on orbitals a or {5 are taken into account. 

This multi-orbital version of the SIAM is also of interest in quantum dot physics, 
where dots with a few orbitals can be prepared and investigated experimentally. 
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One of the interesting aspect of such system is the possibility to observe the orbital 
Kondo effect m. 



3.1 LMA generalization 

In the mean field approximation of the multi-orbital SIAM we also encounter a dou- 
bly degenerate solution, where the two possible Green functions differ only by the 
sign of the impurity magnetic moment. Within the LMA, we introduce for each pair 
of orbital indices a and j3 the two Green functions G"^'^{a))^^ and G"^'^ {co)^^ 
that correspond to the two possible directions of the total magnetic moment on the 
impurity. These Hartree-Fock Green functions depend now parametrically on values 
of local moments on each of the orbitals jia- Next we use the RPA approximation 
to obtain two Green functions G"^'^{a)) and G"^'^{co), which are parametrically 
dependent on the local moments on each orbitals jia- 



3.2 Symmetry restoration and determining the local moment values 

The symmetry restoration in the multi-orbital case is a straightforward generaliza- 
tion of the previous ansatz, i.e. 

Gfico) = \[Gf-\co) + Gf''ico)), (15) 

except that now the symmetrized Green functions G^ {(o) depend explicitly on 
local moments on all of the orbitals, i.e. j/ial- The parameters \iia \ have to be de- 
termined independently. They are found by the minimization of the ground state 
energy of the impurity with respect to both local moment values on orbitals jJ-a and 
particle number on each of the orbitals 

^physical = min folMa^na)- (16) 

As mentioned above, the variational procedure allows us to extend the LMA on the 
multi-orbital cases, where the Luttinger (Fermi liquid) condition for each orbital is 
absent. Also the possibility of non-Fermi liquid solution is naturally included within 
present generalization of the LMA |201, i.e. the variational local moment approach. 
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4 Application to multilevel quantum dots 



A single quantum dot with many atomic-like levels coupled to leads are described 
by a multi-orbital single impurity Anderson model; 

H = //dot + //leads + //dot-leads, 

where //dot is the local impurity part of the SIAM Hamiltonian, //leads corresponds 
to the conduction electron part of the SIAM Hamiltonian, and //dot-leads is equal to 
the hybridization term in SIAM ll24l . 



4.1 V-LMA in quantum dots 

The properties of transport in a quantum dot in equilibrium, i.e. with infinitesimally 
small bias voltage between the leads, are determined by the spectral functions on 
each of the orbitals. Examples of the spectral functions are presented in Fig.|2]for 
the one-orbital case (left panel) and for the two orbital case (right panel). In the two 
orbital case the atomic levels are shifted such that one of the orbitals is at half filling 
(dashed line) and the other is away of half filling (solid line). The Kondo peak in the 
symmetric case is suppressed by the exchange (Hund) interaction (J ^ 0), which 
favors parallel spin orientations. In the asymmetric case the Kondo peak survives 
due to the presence of uncompensated magnetic moment and is shifted toward the 
lower Hubbard band. Further investigation of multilevel quantum dots including 
transport properties will be presented elsewhere l l20l . 




Fig. 2 Spectral functions for one level (left panel) and two level (right panel) quantum dots. Left 
panel: spectral function at half-filling and U = 6 (inset: the ground state energy as a function of 
the absolute value of local moment |^|; the axis starts at the Hartree-Fock value = |/io|)- Right 
panel: orbitally resolved spectral functions in the dot for f/ = 3, 7 = 0.25f/, |ei — e^l = 0. If/, and 
the total filling n,i = 1.95. All curves are for semi-elliptic hybridization function with the width 
W = 20. The Fermi level is at zero. 
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Fig. 3 Spectral functions for two-orbital Hubbard model with different band widths W = 2 (solid 
line) and W = 4 (dashed line) on Bethe lattice with infinite coordination number. Left panel: U = 
1.2, 7 = 0.1(7. Right panel: (7 = 2, / = 0.1(7. The inter-band interaction U' = U -2J is fixed 
preserving 5(7(4) symmetry. In both cases the Fermi level is at zero energy and the bands are half 
filled. 




5 Application to the multi-orbital Hubbard model 

The generalized variational LMA is also applied to solve the multi-orbital Hubbard 
model 

ij a. a 

where the local part is a lattice sum of the terms which are of the same form as the 
atomic part in the SIAM. This model is solved within the DMFT where the self- 
consistency condition relates the local matrix Green functions with the matrix of 
the self-energies lfT4l . In this way the lattice problem is mapped onto the Anderson 
impurity problem which has to be solved for different hybridization functions until 
self-consistency is achieved. In order to solve the Hubbard model within DMFT we 
need to solve the SIAM for arbitrary hybridization functions. The self-consistency 
condition simplifies greatly for the Bethe lattice which is used in this contribution. 



5.1 V-LMA method in DMFT 

In the recent few years the orbital-selective Mott-Hubbard metal-insulator transition 
has been the subject of extensive studies 1221 l25l [T3] l4l l2l . 

Using the V-LMA to obtain the solution of the SIAM in each of the DMFT 
loops the spectral functions for two-orbital Hubbard model at zero temperature were 
found. As an example, Fig.[3]shows the results for the case with different bandwidths 
and non-zero Hund coupling J. Since one of the spectral function is metallic-like 
(finite at o) = 0) and the other is insulating-like (vanishes at © = 0) we conclude 
that the orbital selective MIT occurs in this model system. 

At the end we discuss the V-LMA in perspective to other methods used to solve 
the impurity problem and DMFT equations. The V-LMA belongs to the class of 
approximate, analytical methods like for example the iteration perturbation theory 
(IPT) ITSI . the non-crossing approximation (NCA) lfT2l . or slave-boson theory (SB) 
II29II . and various extenstions of these methods. As we showed here, the V-LMA is a 
conserving approximation, contrary for example to the IPT, and correctly describes 
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high- and low-energy parts of the spectra, recovering the Kondo peak and Luttinger 
pinning. We tested this theory at zero temperature but there is no conceptual obstacle 
why the V-LMA should not work at finite temperatures as well. 

The V-LMA is not numerically exact like the quantum Monte Carlo method lfT4ll . 
the numerical renormalization group (NRG) H], dynamical matrix renormalization 
group (DMRG) |f3T| , or exact diagonalization (ED) (TT] . However, each of the nu- 
merically exact methods suffers from principal obstacles in practical applications, 
in particular when the temperature is too low (QMC) or too high (NRG), or number 
of orbitals increases (NRG, DMRG, ED). 

Therefore we conclude that the V-LMA is a method of choice for solving the 
DMFT equations and can be used as a relatively fast and accurate impurity solver 
The only technical difficulty in the variational LMA is to compute with high accu- 
racy the system energy and to find its minimum. This should be performed with a 
great care. 



Summary 

The generalized variational LMA to the multi-orbital SIAM allows us to efficiently 
solve the problems of correlated electron systems such as multilevel quantum dots 
and the Hubbard model within the DMFT. In particular it is relatively easy to address 
the problems of different band widths and also the removing of the orbital degene- 
racy lfT9l . We experienced that the local moment approach is an efficient method in 
studying these problems, in particular, when the number of the orbitals is larger than 
two. 
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